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to : ABSTRACT 



We present F555W (V), F439W (B), and F336W (U) photometry of 9507 stars 
in the central 2' of the dense, post core collapse cluster M30 (NGC 7099) derived 
from Hubble Space Telescope Wide Field/Planetary Camera 2 images. These 
data are used to study the mix of stellar populations in the central region of the 
cluster. Forty eight blue straggler stars are identified; they are found to be strongly 
concentrated towards the cluster center. The specific frequency of blue stragglers, 
Fbss = N{BSS)/N{V < VuB + 2), is 0.25 ± 0.05 in the inner region of M30 (r < 20"), 
significantly higher than the frequency found in other clusters: -Fbss = 0.05 - 0.15. 
The shape of M30's blue straggler luminosity function resembles the prediction of 
the collisional formation model and is inconsistent with the binary merger model, of 
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Bailyn & Pinsonneault (1995, ApJ, 439, 705). An unusually blue star {B = 18.6, 
B — V = —0.97), possibly a cataclysmic variable based on its color, is found about 
1'.'2 from the crowded cluster center. Bright red giant stars {B < 16.6) appear to be 
depleted by a factor of 2-3 in the inner r < 10" relative to fainter giants, subgiants, 
and main sequence turnoff stars (95% significance). We confirm that there is a radial 
gradient in the color of the overall cluster light, going from B — V ^ 0.82 at r ~ 1' 
to B — V ^ 0.45 in the central 10". The central depletion of the bright red giants is 
responsible for about half of the observed color gradient; the rest of the gradient is 
caused by the relative underabundance of faint red main sequence stars near the cluster 
center (presumably a result of mass segregation). The luminosity function of M30's 
evolved stars does not match the luminosity function shape derived from standard 
stellar evolutionary models: the ratio of the number of bright giants to the number of 
turnoff stars in the cluster is 30% higher than predicted by the model (3.8ct effect), 
roughly independent of red giant brightness over the range My = —2 to +2. 

Subject headings: globular clusters: individual (M30, NGC 7099) — blue stragglers — 
color-magnitude diagrams — stars: evolution — stars: luminosity function 



1. Introduction 

The Hubble Space Telescope {HST) is ideally suited for the study of individual stars in the 
crowded central regions of dense Galactic globular clusters. The projected density of evolved 
stars alone ranges from about 0.3 arcsec"^ in the core of the well studied dense cluster M13 to 
^ 30 arcsec"^ at the center of M15, the Galactic cluster with the highest known central surface 



density ( Cohen et al. 1997 , hereafter Paper VI). Resolving individual stars at these densities is 
beyond the current angular resolution limits of ground-based telescopes, but such observations are 
crucial in order to address several important astrophysical issues. 

Globular clusters are excellent testing grounds for models of stellar evolution because the 
cluster members are coeval and have similar chemical composition (cf. [Bergbusch fc VandenBerg 



1992). Deviations from the natural process of (isolated) stellar evolution caused by stellar 
interactions in these environments of extreme star density can be explored by examining the 
radial dependence of the mix of stellar populations (cf. Piotto et al. 1988| ; Djorgovski et al. 1991 ) 



Additionally, the dense cores of globular clusters are unique laboratories for studying the effects 
of two-body relaxation, equipartition of energy, and binaries on the dynamical evolution of dense 
stellar systems (cf . [Hut et al. 1992 ; pVIeylan fc Heggie 1997| ) . A visible product of stellar collisions 



and mergers of binaries, blue straggler stars (BSSs), are preferentially found in the central regions 
of most globular clusters ( Ferraro et al. 1993| ; |Bailyn Sz Pinsonneault 1995| ). 



This is the seventh in a series of papers describing HST observations of the centers of the 
nearest Galactic globular clusters with |6| > 15°. The main scientific goals of this program are 
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to measure the shape of the density profile in clusters and to understand the nature of evolved 
stellar populations in very dense regions by probing the variation in the mix of stellar types as a 
function of radius (and hence stellar density). Complementary programs targeting main sequence 
(MS) stars in globular clusters are being conducted independently by other groups to explore 
cluster dynamics (cf. |Sosin 1997 ) and models of stellar evolution (cf. piotto et al. 1997] ) . In this 
paper we examine the evolved stellar populations of M30 (NGC 7099) using a set of techniques, 
developed in earlier papers ( Guhathakurta et al. 1992 ; Yanny et al. 199"4a| ; Guhathakurta et al 



1996 , hereafter referred to as Papers I, II, and V, respectively), to: (1) Build empirical point 
spread function (PSF) models using isolated bright stars, allowing for PSF variability across the 
image and using faint stars to reconstruct the saturated cores of brighter stars; (2) Iteratively fit 
the PSF template to the stars on the image; (3) Perform aperture photometry on each star after 
subtracting its neighboring stars with the best-fit PSF template; and (4) Carry out detailed and 
realistic image simulations to assess the effects of crowding on photometric accuracy and sample 
completeness. 

M30 is a prototypical "post-core-collapse" globular cluster based on its rising surface brightness 
profile at radii < 3". Its central surface brightness in the V band is //v'(O) = 15.20 mag arcsec"^ 
( Djorgovski 1993|) , corresponding to a projected density of ~ 5 arcsec"^ in post MS stars alone. 
The cluster has a relatively low metallicity, [Fe/H] = 



2.13 (iDjorgovski fc Piotto 199^ ; |Zinn 



1985 ), a low line-of-sight reddening, Eb-v = 0.05 (|Burstein fc Heiles 1982 ), and is located at a 



distance of about 9.8 kpc ( Reid 1997 ). Based on this recent distance determination, the age of 



M30 is estimated to be about 10 Gyr ( Sandquist et al. 199^ ). It has long been noted that M30 



displays a central bluing trend like other post-core-collapse clusters (Williams 



ers (Wihiams & Bahcah 1979|; 


Piotto et al. 1988 




Burgarella & 



Buat 19961). 



In this paper we use HST Wide Field/Planetary Camera 2 (WFPC2) data to investigate the 
origin of M30's nuclear color gradient and to study its abundant population of BSSs; Yanny et al. 



(1994b, hereafter Paper IV) used the same data set to examine the inner 20" of M30. This paper 
is organized as follows: Sec. 2 contains a description of the data and reduction techniques; Sec. 3 
contains results from our study of M30's evolved stellar populations, with particular emphasis 
on BSSs, radial population gradients, and the stellar luminosity function (LF); Sec. 4 contains a 
summary of the main points of the paper. 



2. Data 
2.1. Observations 

A set of eight WFPC2 images of the center of M30 (NGC 7099) were obtained on 1994 
March 31 (see Paper IV for a discussion of the data set): 2 x 100 s with the F336W filter, 2 x 40 s 
with F439W, and 4 x 4 s with F555W. These filters are roughly similar to the Johnson U , B, and 
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V bandpasses, respectively. All eight exposures were obtained at the same telescope pointing and 
roll angle. 

Each WFPC2 image consists of four 800 x 800 Charge Coupled Device (CCD) frames: a 
Planetary Camera CCD (PCI) with a scale of 0'.'0455 pixel"^ and three Wide Field Camera 
(WF2-WF4) CCDs with a scale of 0'.'0996 pixel'^ The usable field of view is 34" x 34" for PCI 
and about 76" x 76" for each WF CCD; the total image area is about 5.1 arcmin^. An image of 
the cluster is shown in Fig. |^ (Plate XXX). The telescope pointing was chosen so that the center 
of the cluster was positioned near the center of the PCI CCD, shown in the lower right quadrant 
of the greyscale WFPC2 mosaic image. The WF2~WF4 images are oriented counterclockwise 
starting from the upper right quadrant of the mosaic. A detailed description of the instrumental 
parameters and in-orbit characteristics of WFPC2 is given by [Burrows (1995) and by Trauger 
et al. (1994)|. 



The exposure times have been chosen so that the brightest stars {V ~ 11.8, B ~ 13.3, 
U ~ 14.5) in the image are not saturated by more than a factor of 4-5 in each individual 
exposure. In general, no more than 10 pixels are affected per star, and it is possible to recover 
their magnitudes and positions using information in the unsaturated wings of the PSF. 



2.2. Data Processing / Stellar Photometry Technique 

We have developed and tested procedures for deriving stellar photometry from HST images 
of dense globular clusters. The reader is referred to Papers I and II for a detailed description. 
Even though the procedures were originally designed for pre-repair HST data with an aberrated 
PSF, several of the steps are relevant for analysis of crowded star fields and we apply them here. 



The following steps for processing WFPC2 data are summarized below (see Guhathakurta et al 



19961 ): cosmic ray removal, image coaddition, star finding, PSF building, PSF fitting/subtraction. 



and aperture photometry. 

Standard bias and flat field calibrations have been applied to the data as part of the 
preprocessing pipeline at the Space Telescope Science Institute. Pixels affected by cosmic rays 
are identified and masked using an algorithm that compares each image to the other image(s) in 
that band; unmasked pixels in all images in a given band are averaged to obtain the final F336W, 
F439W, and F555W images. 

Cross-correlating pairs of images in different bands shows that the stellar positions in the 
F555W and F439W band PCI images are offset from the corresponding stellar positions in 
the F336W band PCI image by 1.2 pixel (0'.'055) and 0.7 pixel (0'.'032), respectively. For the 
WF2-WF4 CCDs, the offset is in the range of 0.34-0.59 pixel (0'.'034-0'.'059) between the F555W 
and F336W images and in the range 0.15-0.32 pixel (0'.'015-0'.'032) between the F439W and 
F336W images. The combined F336W, F439W, and F555W images are aligned by fractional 
pixel interpolation and summed to produce a deep image. Interpolation results in degradation 
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Hubble Space Telescope WFPC2 Image of the Central Region of IV130 




Fig. 1. — A F-band WFPC2 image of M30 (negative greyscale) and a "true color" U + B + V 
composite of the cluster's central 15" x 15" (color inset), along with the image scales and orientation. 
The cluster centroid is indicated by the green '+'. The colored circles indicate: 29 blue stragglers 
(blue); a very blue star near the cluster center (mauve); and the astrometric reference star near the 
top of the color inset (red). 
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of the PSF by an amount that varies from star to star due to the undersamphng of the PC and 
particularly the WF CCD data; however, the interpolated image is only used to find stars (see 
below) and not to derive stellar photometry. 

A preliminary list of stars is obtained by applying the standard peak finding algorithm find 
of the DAOPHOT package ( ^tetson 198^ ) to the summed F336W+F439W+F555W image. Bright 
relatively isolated stars are used to build an empirical, quadratically variable daophot ii PSF 



template ( Stetson 1992 ) for each of the four bands and for each of the four CCD frames (PCI, 
WF2-WF4). This PSF template is fit to and subtracted from the stars on the image using 
DAOPHOT ii's ALLSTAR routine. After inspecting the residual images (original minus best-fit 
template), the star list is edited to remove spurious objects, such as PSF artifacts around bright 
stars and hot pixels (the slight offset between stellar positions in F336W and F555W images 
facilitates the identification of fixed pattern artifacts), and to add stars fainter than the find 
threshold. The process of PSF building, fitting, and subtraction is iterated a few times; in each 
instance, the results of the previous iteration are used to remove neighbors of the bright stars 
used for PSF building. Finally, aperture photometry is obtained for each star, after its neighbors 
have been removed using the best-fit PSF template. The photometry lists in the three bandpasses 
are merged by positional matching. The final list of stars identified and matched in each of 
the three filters consists of 3114 stars from the PCI image (0.32 arcmin^), and 2055, 1495, and 
2843 stars from WF2-WF4, respectively (each WF CCD covers 1.60 arcmin^). 

Tests based on bright, relatively isolated stars show that the PSF is stable across the WFPC2 
field of view (see also Papers IV and V): the correction term for converting r ~ O'.'l aperture 
magnitudes to r = 1" "total" magnitudes varies by a few percent or less across each CCD image. 
Furthermore, this slight variation is similar in all three bands and the telescope pointing is 
identical for all exposures in all bands so the effect on measurement of stellar colors is negligible. 

Differences in the zeropoints of the magnitude scale of the individual CCDs are determined by 
cross-correlating the "positions" of bright {V ^ 18) stars in the color-magnitude diagram (CMD) 
between pairs of CCD frames. The tightness of the CMD features (Fig. |^, and the fact that 
several hundred stars are used to define the peak of the cross-correlation function, makes it possible 
to determine color and magnitude offsets to an accuracy of ^ 0.01 mag. The inter-CCD magnitude 
zeropoint offsets for WF2-WF4 with respect to PI are [+0.25, +0.25, +0.25] in F555W, [+0.23, 
+0.22, +0.20] in F439W, and [+0.43, +0.39, +0.41] in F336W. The zeropoint difference between 
the WF CCDs and PCI results primarily from the different aperture sizes used to derive stellar 
photometry: rwF = 0'.'16 vs rpci = 0'.'09. The zeropoint difference is largest in F336W because its 
PSF core is significantly broader than in F439W and F555W. 

The cross-correlation technique used for inter-CCD zeropoint matching, as well as the 
procedure used to match the WFPC2 photometry to ground-based photometry (see below), relies 
on the validity of an implicit assumption: That there are no intrinsic spatial gradients in the color 
of the individual evolved stellar populations in M30, particularly in the B — V and U — V colors 
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of faint red giant branch (RGB) stars which dominate by number the V ^ 18 sample used for 
matching. Radial gradients in RGB color have never been convincingly demonstrated to exist in 
globular clusters; it is common to find changes in the mix of stellar types, but not in the color of 
a given stellar type. Moreover, it is reassuring that there is a perfectly plausible explanation for 
the observed inter-CCD zeropoint differences in the M30 data set. 



A charge transfer efficiency problem during WFPC2 readout (Burrows 1995; Holtzman] 



et al. 1995b ) causes objects at high row number to appear systematically fainter than if they had 
been at low row number. For exposures taken at — 77°C, a linear ramp correction ranging from 
no correction at row 1 to —0.1 mag at row 800 (original CCD coordinates) works well to first 
order. We apply such a linear ramp correction to the stellar photometry; it results in a noticeable 
tightening of the horizontal branch (HB) and subgiant branch/main sequence turnoff (MSTO), 
the only features that are more or less horizontal in the CMD. Note that color measurements are 
unaffected by the charge transfer efficiency problem since the telescope pointing was identical for 
exposures in all three bands. 

The zeropoints of the WFPC2 instrumental F336W, F439W, and F555W magnitudes are 
adjusted to roughly match (~ 0.1 mag) ground-based photometry of stars in M30 in the Johnson 
U, B, and V bands ( Dickens 1972| ). There is not sufficient overlap between this ground-based 



UBV data set (or any other published ones) and the HST data set to allow an extensive 
star-by-star comparison, the cluster's central region being too crowded at arcsecond resolution. 
Instead, various fiducial points in the CMD — HB brightness, and the RGB color at the level of 
the HB and at the RGB base — and a variety of filter combinations [{B — V^ V) and (U — B, B)] 
are used to determine the zeropoints in all three bands. These zeropoints are in good agreement 
with the published WFPC2 photometric zeropoints (Burrows 1995). Following the prescription 
of Holtzman et al. (1995a)| , the zeropoint-adjusted magnitudes on the WFPC2 bandpass system, 



m336, ?Ti439, and msss, are converted to the Johnson UBV system: 

U = mssQ - 0.240 (msse - m^^^) + 0.048(^336 -171555)'^ 

B = 771439 + 0.003(777439 -777555) - 0.088(777439 -r77555)^ 

V = 777555 - 0.060(777439 -777555) + 0.033 (7r7439 -"7555)^ (1) 



The B and V magnitude conversions are accurate, but the converted "f7" magnitude may 
be significantly different from true Johnson U because of the severe red leak in the F336W 
bandpass (e.g., see Fig. |6|) which is not adequately accounted for by quadratic terms in the color 
transformation equation. 
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2.3. Astrometry 

Astrometric calibration has been carried out using the metric task (Gilmozzi et al. 1994) 
within iraf's stsdas package. This task uses HST pointing information based on the Guide 
Star Catalog to derive absolute astrometry, and corrects for mechanical geometric distortion and 
spherical aberration that vary from CCD to CCD, to place all the stars on a global coordinate 
system. A relatively isolated bright star, ID# 3611 in Table ^{V = 15.51, B = 16.14, U = 16.00), 
is adopted as the astrometric reference point: 



aREF(J2000) = 2r40""22';402, dREF(J2000) = -23°10'42'.'42 . (2) 

The astrometric reference star (red circle in Fig. |l]) is located 2'.'3 E and 7'.'6 N of the cluster 
centroid (green cross in Fig. |l]). 



The number-weighted centroid of evolved stars in M30, as determined in Paper IV , has 
coordinates of: 



ao(J2000) = 2r40™22!24, 5o(J2000) = -23°10'50'.'0 . (3) 

to a Id accuracy of 6r = 0'.'2. We adopt this as the cluster center in the rest of the analysis 
presented in this paper. 



2.4. Photometric Accuracy and Completeness 

Using artificial star tests and realistic simulations of HST WFPC2 images of M30 and M15 
carried out in Papers IV and V, respectively, we estimate the degree of completeness and level of 
internal photometric accuracy in the data set analyzed in this paper. These error estimates, as 
well as the other error estimates discussed below, refer to the internal accuracy in the measurement 
of relative instrumental magnitudes; there may be additional systematic errors in the conversion 
to Johnson magnitudes (^ 0.1 mag in B and V, and significantly larger in U due to the F336W 
red leak). The ratio of exposure times for M30 and M15 has been chosen to compensate for 
the 0.76 mag difference in their distance moduli ( Pjorgovski 199^ ) so that the ratio of signal to 
(read+Poisson) noise is the same for stars of a given absolute magnitude. The surface density of 
post MS stars is about 10 times higher in M15 than in M30 (at comparable angular distances from 
the cluster center), while the shape of the LF of evolved stars is similar in the two clusters. Thus, 
the effect of crowding is slightly less severe in the M30 data set than in the M15 data set — i.e., the 
degree of completeness is higher and the photometry more accurate in M30 at a given absolute 
magnitude and angular distance from the cluster center. 

Our best estimate for the la photometric error in V in the dense central 10" of M30 
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is ^ 0.05 mag for bright stars (V ^ 17) and about 0.10-0.15 mag for stars at the brightness level 
of the MSTO (V ^ 18.6). In this region of M30, effects of blending are expected to start becoming 
noticeable around the MSTO. The sample should be complete down to the MSTO even in the 
crowded central regions of the cluster. The simulations also indicate that, within r < 5", the level 
of completeness starts falling off for stars just below the MSTO {V ^ 19), dropping to roughly 
50% for stars 2 mag below the MSTO. Incompleteness is expected to set in at progressively fainter 
magnitudes for samples at larger radii. 

We carry out a consistency check of these completeness estimates by studying the observed 
LF of M30 stars as a function of radius. The stars detected in the full WFPC2 mosaic are 
divided into eight radial bins around the cluster center: (1) r < 5'.'00, (2) 5'.'00 < r < 9'.'80, 
(3) 9'.'80 < r < 15'.'41, (4) 15'.'41 < r < 23'.'2, (5) 23'.'2 < r < 35'.'8, (6) 35'.'8 < r < 51", 
(7) 51" < r < 71", and (8) 71" < r < 130". The boundaries of these radial bins are chosen so 
that each bin contains roughly the same number of stars with msss < 18.6 (~ Vmsto); the data 
set should be complete at all radii for this bright subsample of stars so it provides a convenient 
means of matching the relative normalization of the LFs in the various radial bins. The shape of 
the LF of post MS stars is expected to be the same at all radii. This is because stars brighter 



than My = M^^TO ^ .^3 5 within a very narrow mass range (AM ~ 0.03 Mq — [Bergbusch &:| 



VandenBerg 1992 ) so that the effects of mass segregation are negligible (cf. pryor et al. 1986 



Bolte 1989). 



The F555W-band stellar LFs are shown in Fig. ||. Note, these LFs are based on stars detected 
on the F555W image, without requiring that the stars also be detected in F439W or F336W; 
hence Johnson V magnitudes are not available for all these stars. The bold line with dots shows 
the "most complete" combined LF based on the average of all radial bins for V < 18 (to minimize 
Poisson error) and on the average of the two outermost annuli (r > 51") for stars with V > 18 
(since the degree of completeness is highest in this sparse outer region). The LF in the innermost 
bin (r ^ 5") is identical to the combined LF down to the MSTO, but peels away at fainter 
magnitudes as expected, dropping sharply beyond V ~ 20.5 as a result of incompleteness. The 
V magnitude at which incompleteness sets in appears to increase roughly monotonically with 
increasing annular radius in keeping with the trend seen in the simulations described above. 

The F336W-band stellar LFs shown in Fig. |^ have a qualitatively similar behavior to those in 
the F555W band. A notable difference between the two bands is that RGB stars span a narrower 
range of brightness in U (3.5 mag) than in y (5-6 mag), a natural consequence of the fact that 
bright RGB stars tend to be redder than faint ones. Thus, the brightness contrast between faint 
MSTO stars and the brightest stars (at the tip of the RGB) is smaller at the shorter wavelength, 
so MSTO stars are slightly less affected by crowding in U than in V. This effect is partly offset 
by the fact that the signal-to-noise ratio for faint MS stars is slightly higher in V than in U: the 
longer exposure time is not enough to compensate for the lower instrumental efficiency and fainter 
intrinsic stellar brightness in U. Another difference between the LFs in the two bands is that the 
slope of the MS portion of the LF is shallower in U than in V; this is because faint MS stars tend 
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Fig. 2. — The observed stellar luminosity function of M30 in the F555W band (roughly Johnson 
V) in eight radial bins. The rollover beyond F555W ~ 21 is caused by incompleteness, which sets 
in at brighter magnitudes at smaller radii because of increased crowding. The bold solid curve 
with dots represents the "most complete" combined LF based on the average all radial bins for 
F555W < 18 (to minimize Poisson error) and on the average of the outermost two bins for stars 
with F555W > 18 (since the degree of completeness is highest in this sparse outer region). The 
main sequence turnoff (MSTO) at F555W = 18.6 is marked; the bump at F555W ~ 15-16 is 
caused by stars in the HB. 
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Fig. 3. — Same as Fig. |2| for the F336W band (roughly Johnson U). The rollover beyond 
F336W ~ 22 is caused by incompleteness. As in Fig. ^, the "most complete" combined LF (bold 
solid line with dots) is based on the average of all bins for F336W < 18 and on the average of the 
outermost two bins for F336W > 18. The main sequence turnoff (MSTO) at F336W = 18.5 is 
marked; the bump at F336W = 15.5 is produced by HB stars. The LF slope in the RGB region is 
steeper in the U band than in V, while the MS slope is shallower, as expected from the orientation 
of the RGB and MS tracks in the CMD. 
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to be redder than bright ones. We return to a detailed discussion of the stehar LF in Sec. 3.5. 



The widths of prominent features in different CMDs (Figs. ^ ^, and |6|) — the RGB at 
= 15 - 17 and the MS just below the MSTO {B = 19 - 20) — are used to estimate the overall 
error in stellar photometry, under the assumption that the features have zero intrinsic width. 
These error estimates represent averages over the entire area of the WFPC2 mosaic. For the filter 
combinations, U — B, B — V, and U — V, the rms scatter in RGB color about empirically defined 
fiducial lines is 0.04, 0.02, and 0.04 mag, respectively, and the rms scatter in MS color is 0.08, 0.08, 
and 0.08 mag, respectively. These translate to photometric errors (in magnitudes) of ajj = 0.03, 
as = 0.01, and ay = 0.02 for the bright stars comprising the RGB, and au = 0.05, as = 0.06, 
and ay = 0.05 for faint MS stars. 



3. Stellar Populations 

Our analysis of M30's stellar populations is based on photometry derived from the HST 
WFPC2 images of the central region of the cluster. Table |l| lists the relative positions, Aq(J2000) 
and A(5(J2000), projected radial distance from the cluster centroid, and Johnson UBV magnitudes 
and colors of the 40 brightest stars located within r < 10" of M30's center. Stellar positions are on 
the equinox J2000 coordinate system and are measured in arcseconds relative to the astrometric 



reference star, ID# 3611, whose coordinates are given in Eq. (2) in Sec. 2.3. The coordinates of 
the (number- weighted) cluster centroid are given in Eq. (3). The complete electronic version of 
Table |l] is available in the electronic edition of the Journal by link to a permanent database and 
via anonymous ftp (see Appendix); it contains 9940 stars from the full M30 WFPC2 data set. 
The stellar ID^s run sequentially through the full table in order of increasing right ascension. 



3.1. Color-Magnitude Diagrams 

Figure § shows a {B — V, B) CMD of all stars detected {B <^ 23) in the 5.1 arcmin^ area of 
the WFPC2 mosaic image of M30. The 9507 stars matched in the U, B, and V filters are marked 
as squares. An additional 433 stars that are detected and position-matched only in B and V are 
marked as dots (most are fainter than V = 21). This CMD is used to assign a stellar type to 
each star: bright RGB and asymptotic giant branch (AGB) stars, faint RGB and subgiant stars, 
HB stars, and blue straggler stars (BSSs). The HB is entirely on the blue side of the RR Lyrae 
instability strip, with its S-band brightness becoming fainter towards bluer B — V colors. The 
bright part of M30's RGB and AGB are relatively blue (relative to the average Galactic globular 
cluster), in keeping with its low metallicity: [Fe/H] = —2.13 ( [Djorgovski 199^ ). The MS is 
reasonably well defined up to ~ 2 mag below the MSTO. There is a prominent, well-populated 
blue straggler sequence in the diagram (Sec. ^^ ). The definition of BSS is somewhat arbitrary, as 
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Table 1. Stellar photometry and astrometry'' for M30'^ (40 brightest stars within r < 10"). 
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^Positions relative to star ID# 3611; see Sec 2.3 for details. 

^The complete table is available in the electronic edition of the Journal by link to a 
permanent database and via anonymous ftp (see Appendix for details). 
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Fig. 4. — A B vs B — V color-magnitude diagram of all 9507 stars in the WFPC2 image of M30 
that are detected and position- matched in the F555W, F439W, and F336W bands bands (squares). 
The 433 stars matched only in F555W and F439W are plotted as dots. The dashed lines indicated 
the boundaries that have been used to define various stellar types. There are 48 stars in the region 
labeled BSS. The large diamond symbol marks an unusually blue star located close to the cluster 
center. The five encircled symbols represent stars which are significantly redder than the red giant 
branch in all three color-magnitude plots (see Figs. |5| and these may be Galactic field dwarfs in 
the foreground. 
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is the boundary between faint and bright RGB stars. The bright vs faint RGB boundary has been 
drawn at B = 16.6, fainter than the level of the HB (the customary demarcation line), in order to 
ensure that the bright RGB/AGB sample is large enough for a study of radial gradients. Note, 
a similar CMD (based on PCI data alone) is presented in Paper IV| . Even though the x-axis of 



Fig. 4 of Paper IV is mistakenly labeled as "C/ — y , the quantity that is plotted in that figure is 



B — V color with an arbitrary offset. 

Using different symbols for different stellar types [as defined by the {B — V, B) CMD], we 
plot the same 9507 stars in the {U - V, B) and {U-B, U) CMDs in Figs. | and ||, respectively. 
The long U — V color baseline is sensitive to subtle differences in stellar temperature (despite the 
red leak in the HST F336W filter), resulting in a clear separation among the HB, bright RGB, 
and AGB. The MS and BSS sequence are mostly vertical in the {U — V, B) projection making it 
difficult to distinguish between "true" BSSs and artifacts caused by the blending of MSTO stars. 
The HB stars all have about the same U brightness, the bolometric correction to the U magnitude 
being nearly independent of temperature for these hot stars, so that the HB feature is nearly 
horizontal in the {U — B, U) CMD. This projection highlights the seriousness of the red leak in 
the F336W filter which causes the bright RGB and AGB to shift to the left so that they intersect 
the HB. The red leak also causes the faint RGB and MS stars to be mapped all the way across to 
the left (blue) side of the BSS sequence. Even with the red leak however, stars of a given type, as 
defined in the {B — V, B) CMD, tend to remain grouped together in the (U — B, U) projection. 

Note the very blue star in the {B-V, B) CMD: B - V -0.97, B ~ 18.6 (large diamond 
symbol in Fig. |; star ID# 2931 in Table |). This star may possibly be a cataclysmic variable star 
based on its location in the CMD, although we are unable to comment on its variability from our 
data set alone. It is located 1".2 SSW of the cluster center, well outside the r = 0'.'6 (formal) 3(T 
error circle for the cluster centroid, offset by Aa = — 2'.'88 and = — 8'.'59 from the reference star 



(see Sec. 2.3 for details on astrometry); this unusually blue star is marked by a mauve circle in 
Fig. 1^. The apparent proximity of the object to the cluster center makes it particularly interesting, 
but the extreme crowding in this region makes spectroscopic followup a difficult prospect (the star 
has a Ay = 5 mag brighter RGB neighbor only 0'.'4 away). 

There are five relatively bright {V ~ 15.5 - 18) objects in the CMDs which are significantly 
redder \\iU — B and B — V (by up to +1 mag) than typical RGB, subgiant, and MSTO stars of 
comparable apparent brightness in M30; these objects have stellar ID#s of 890, 4735, 9200, 9350, 
and 9508 in Table || and are indicated by encircled symbols in Figs. 1^-^. One cannot rule out 
the possibility that these are cluster stars with abnormal spectral characteristics. However, the 
fact that these five objects are not concentrated towards M30's center casts serious doubt on their 
cluster membership: four out of the five lie beyond r = 1', whereas only about 15% of cluster 
members within the WFPC2 image are expected to lie beyond this radius. No radial velocity 
measurements are available for these five red objects (K. Gebhardt, private communication); 
M30's large negative radial velocity, fMso ~ —187 km s~^ ( [Zaggia et al. 199"2[ ), should make it easy 
in principle to distinguish field stars from cluster members. On the basis of a realistic model of the 
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Fig. 5. — Same as Fig. ^ except B magnitude is plotted against U — V color for the 9507 stars 
matched in all three bands. Different symbols are used to plot the different stellar types, defined 
according to the {B — V, B) CMD (Fig. ^). This projection distinguishes clearly between the 
HB, bright RGB, and AGB with the help of the wide U — V baseline which is sensitive to subtle 
differences in stellar temperature. 
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Fig. 6. — Same as Fig. |5| except U magnitude is plotted against U — B color. The HB is nearly 
horizontal in this projection. The red leak in the WFPC2 F336W filter causes the RGB and MS to 
be artificially shifted to the left, so that they overlap with the HB and BSS sequence, respectively. 
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Galaxy, Ratnatunga fc Bahcall (1985)| predict a surface density of field star interlopers of about 
3 arcmin"^ to V ^ 23, or a total of 15 interlopers within the area of the WFPC2 mosaic. Most of 
these field stars are expected to be faint: less than four field stars are expected to be brighter than 
V ^ 21 and less than one star (0.37) is expected to be brighter than V = 19. The fact that even 
field stars tend to cluster together implies that the probability of finding five bright field stars 
could be higher than the Poisson estimate. If these five red objects in the the M30 WFPC2 data 
set are foreground MS stars of similar metallicity, there is a 3 mag spread in distance modulus 
(factor of 4 in distance) amongst them. These objects are even less likely to be background 
galaxies. Field stars are expected to outnumber field galaxies for V < 20, even at high Galactic 
latitude. Moreover, our PSF fitting procedure is designed to identify and eliminate extended 
objects (galaxies); no such objects are found in the M30 data set. 



3.2. Blue Stragglers 

3.2.1. Identification 

The {B — V, B) CMD of M30 presented in Fig. |^ shows a large collection of BSS candidates 
in the region between the MSTO and the extended blue tail of the HB. This sequence of BSS 
candidates lies roughly where the continuation of the MS would be, if it were to be extrapolated 
beyond (i.e., brighter and bluer than) the MSTO. These stars are, for the most part, distinct from 
HB and MSTO stars in all three color-magnitude projections (see Figs. §-P). 

Table § lists the positions, brightnesses, and colors of the 48 BSS candidates outlined by the 
dashed lines in Fig. ^. All positions are on the equinox J2000 coordinate system and are measured 



with respect to the reference star as described in Sec. 2^. The stars in Table g are arranged 
in order of increasing projected distance from the cluster center. The stellar ID^^s for the BSS 
candidates are the same as those in the complete (electronic) version of Table ||. 

A "f " symbol is used to mark marginal candidates in Table |2|, those for which the BSS 
designation may be in some doubt. Results from earlier image simulations (Papers IV and V) 
indicate that blending and/or errors in the photometry of subgiant/MSTO stars can contaminate 
the BSS portion of the CMD. We have been very conservative in designating "marginal" BSS 
candidates: these include all the ones fainter than i? = 18 (since this region may be populated 
by blend artifacts), those that are close to the edge of the selection region, and outliers that are 
separated from the upward extrapolation of the MS track. It is possible that some fraction of 
the marginal BSS candidates in M30 are actually subgiant/MSTO stars that have been scattered 
away from their fiducial locations in the CMD due to measurement error. Contamination is likely 
to be most severe at the base of the BSS sequence where it meets the subgiant branch and MSTO 
[B f» 18.5, B — V 0.25), so this region has deliberately been excluded from the BSS defining 
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Table 2. Blue straggler candidates in M30. 
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+ Marginal blue straggler candidates (see Sec. 3.2.1 for details). 



- 20 - 



box from the outset (Fig. Our BSS sample is likely to be incomplete at the faint end (near the 
MSTO), as the defining box is liable to exclude BSS if their luminosity -Less ^ 1-4Lmsto- The 
exact balance between the competing effects of BSS sample contamination and incompleteness 
depends on the details of the somewhat arbitrary BSS defining criteria used. The M30 BSS 
population is unlikely to be contaminated by field star interlopers; less than one field star is 
expected in the region of the {B — V, B) CMD brighter and bluer than the MSTO based on the 
Bahcall-Soneira model of the Galaxy ( Bahcall &: Soneira 19801 ; Ratnatunga & Bahcall 1985). 



3.2.2. Specific Frequency 



The specific frequency of BSSs in M30 is calculated according to the definition of [Bolte et al 



(1993)| : Fbss = N{BSS)/N{V < I^hb + 2). Although M30's HB stars do not ah have the same 
V brightness (Fig. the HB has a well defined brightness of ^ ~ 15.1 at the location of the 
RR Lyrae instability region {B — V^ 0.3) before dropping to fainter V magnitudes at the blue 
end. There are 332 RGB/AGB/HB/bright subgiant stars in the 5.1 arcmin^ WFPC2 mosaic image 
(r ^ 130") that fulfill the criterion for the normalizing population and 48 BSSs — an overall specific 
frequency of -Fbss = 0.14 it 0.02 {la Poisson error). This statistic is computed as a function of 
radius and the results are listed in Table |3[ The numbers listed in Table |3| are slightly different 



from those quoted in Paper IV; this is due in part to the inclusion of data from the WF CCDs and 



in part to a slightly revised defining criterion for BSS candidates. 



Table 3. Blue straggler and subgiant /turnoff specific frequency vs radius in M30. 
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The BSS specific frequency found in the inner 20" of M30, Fbss(^ < 20") = 0.25 ± 0.05, is 
higher than the specific frequency found in other well studied (e.g., at HST resolution) globular 
clusters, including ones whose density is comparable to or greater than that of M30, such as M15 
(cf. ISosin k King 1995| ; Paper VI). In fact, the overall BSS frequency in the M30 WFPC2 data 
set, -Fbss = 0.14, falls at the high end of the typical factor of 3 range seen in Galactic globular 
clusters. 
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If the marginal BSS candidates in M30 are excluded from the analysis, the remaining (secure) 
BSS candidates have a specific frequency of -Frss = 0.08 in the overall data set and -Fbss = 0.13 in 
the central r < 20". These values should be treated as lower bounds to the "true" BSS frequency 
in M30 — a fair fraction of the marginal candidates (possibly the majority) are likely to be real 
BSSs, and the defining box is liable to have missed faint BSSs (Sec. \3. 2. i| ) . For the purpose 
of comparing M30's BSS frequency to that of other clusters, it is best to consider Fbss values 
based on the full set of BSS candidates (secure+marginal) because: (1) estimates of Fbss in other 
clusters tend to include marginal candidates, and (2) the BSS defining criteria used in this paper 
(Fig. ^) are similar to those used in other studies. 



3.2.3. Central Concentration 
The value of Fbss ™- M30 displays a dramatic increase towards the cluster center (see Table ^ 



and Sec. 3.3.1 ). The BSS frequency in the 20" < r < 130" region of the WFPC2 mosaic image 
is 0.06 lb 0.02, a factor of 4 lower than in the inner r < 20" (Fbss = 0.25), but in keeping with the 
value seen in other clusters (Paper VI). The concentration of the M30 BSSs towards the cluster 
center is consistent with the idea that they are more massive than individual RGB/subgiant stars 
(either merger remnants or members of close binary systems), and have consequently sunk deeper 



into the cluster potential well (King 1991). 



Restricting the M30 BSS sample to only secure candidates (i.e., those not marked by a 
"I" in Table ^, the specific frequency is seen to increase by a factor of 3 from Fbss = 0.04 in 
the outer 20" < r < 130" region to Fbss = 0.13 in the central r < 20" region. The sample 
of secure BSS candidates, being well removed from the subgiant/MSTO region of the CMD, is 
uncontaminated by blends. This is firm proof that M30's BSSs are indeed physically concentrated 
towards the cluster center; the observed concentration is not an artifact of increased photometric 
scatter /blending in the crowded central region of the cluster. In fact, the data in Tables |2| and 
^ suggest that the majority of even the marginal candidates in M30 are faint (but real) BSSs 
which are centrally concentrated in a similar manner to their brighter (and hence more secure) 
counterparts. 

We investigate whether the high BSS fraction in M30's inner region is produced by anomalies 
in the normalizing population of RGB/ AGB/HB/bright subgiant stars (V < Vhb + 2) by 
comparing the radial distribution of this normalizing population to that of faint subgiant/turnoff 
(SGTO) stars in the cluster. The SGTO class is defined by the criterion: Vhb + 2 < V < FmstO) 
where Vmsto = 18.6. The fraction, Fsgto = N{V < Vhb + 2)/A^(Vhb + 2 <V < Fmsto), is 
consistent with being constant throughout the central r < 2' of M30 (Table |3|); the apparent 
5%-10% decrease in Fsgto in the inner 20" compared to the global value is not statistically 
significant (^ la effect), and this has a negligible effect on the measurement of Fbss- Furthermore, 
a study of the LF of M30's evolved stellar population shows that the abundance of RGB/AGB 
stars relative to fainter stars (faint subgiant, MSTO, and MS stars) is not anomalously low; rather, 
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their abundance is a little higher than normal (see Sec. |3.5[). 



3.2.4- Blue Straggler Luminosity Function 

Theoretical studies of the formation of BSSs indicate that the BSS LF may be a good 
discriminant between different formation scenarios ( Bailyn &: Pinsonneault 1995| ). For example, 



stellar collisions produce a tail of very luminous BSSs while mergers of primordial binaries do 
not. Bailyn & Pinsonneault adopt a top heavy input mass function for stellar collisions since 
these are important in the dense central regions of clusters where massive stars are expected to be 
concentrated as a result of mass segregation; input masses for the primoridal binary merger case 
are likely to be drawn from a more normal (e.g., ^alpeter 1955| ; x = +1.35) mass function. The 
authors also consider two extreme options for the compositional structure of a newborn BSS that 
ought to bracket the realistic range of possibilities — they dub these 'collisional' chemistry, which 
corresponds to complete mixing yielding a chemically homogenous product, and 'binary merger' 
chemistry, which corresponds to no mixing. 

Figure ^ shows the observed BSS LF in M30 (bold dots with Poisson error bars) compared 
to various Bailyn & Pinsonneault model BSS LFs for metal poor stars: Z = 6.25 x 10~^ or 
[Fe/H] = —1.7. The BSS luminosities are estimated from their U brightnesses, with a +0.25 mag 
correction to the observed f/MSTO = 18.5 to account for red leak and an additional —0.1 mag 
differential bolometric correction. We expect our sample of M30 BSSs to be incomplete for 
-^BSsZ-^^MSTO < 1-4 or log(-LBss/-^MSTo) < 0.15 (indicated by vertical dot-dashed line in Fig. |^); 
this region is too close to the MSTO and our BSS selection boundary has been designed to avoid 
contamination from MSTO and subgiant stars. The model LFs and the M30 data are normalized 
to unity over the "complete" range, log(L/LMSTo) ^ 0.15. 

The shape of the LF of BSSs in M30 bears a general resemblance to the shape of the collisional 
model LFs (x = —2). The peak in the M30 LF at log(LBss/-^^MSTo) ~ 0.5 and the suggestion of 
a tail to higher BSS luminosities are not reproduced in the primordial binary merger model LF: 
An X = +2 input mass function binary merger model is shown in Fig. ^; in fact, all merger models 
with X 1 have a qualitatively similar shape to the x = +2 case ( Bailyn &: Pinsonneault 1995| ). 



The intermediate luminosity peak in M30's BSS LF is also apparent from the {B — V, B) CMD 
in Fig. ^ note the concentration of BSSs and the hint of a "gap" between the BSS sequence and 
MSTO. While BSSs are more abundant in M30's central region than in any of the other clusters 
studied to date, the total number of stragglers found in the WFPC2 image of M30 is too small to 
draw any firm conclusions about the shape of their LF. It is fair to say though that the data are 
strongly suggestive of a collisional origin for the BSSs. 
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Fig. 7. — Luminosity function of the blue straggler stars in M30 (bold dots with Poisson error 
bars). The BSS sample is likely to be incomplete fainter than log(L/LMSTo) = 0.15 (dot-dashed 
vertical line) due to our selection criterion. Three theoretical LF models (Bailyn &; Pinsonneault 
1995) are plotted for comparison, reflecting a variety of BSS formation scenarios: x = —2 vs x = +2 
input mass function slopes (appropriate for stellar collisions and binary mergers, respectively), and 
complete mixing vs no mixing for the chemical composition of the newborn BSS. The models are 
based on a metal abundance of Z = 6.25 x 10~^ ([Fe/H] = —1.7) and, like the data, are normalized 
to unity over the range log(L/L]v[STo) > 0.15. 
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3.3. Population Gradients 

In order to study radial variations in the mix of M30's stellar populations, the stars detected 



in the full WFPC2 mosaic are divided into eight annuli around the cluster center, as in Sec. 2.4: 
(1) r < 5'.'00, (2) 5'.'00 < r < 9'.'80, (3) 9'.'80 < r < 15'.'41, (4) 15'.'41 < r < 23'.'2, (5) 23'.'2 < r < 35'.'8 
(6) 35'.'8 < r < 51", (7) 51" < r < 71", and (8) 71" < r < 130". The limiting radii of the 
eight radial bins have been chosen so that each bin contains approximately the same number of 
stars with V < 18.6, a limiting magnitude for which the samples at all radii are expected to be 
complete. 



3.3.1. Comparing the Radial Distribution of Various Stellar Types 

In this section, we compare the radial density distributions of the various types of evolved 
stars found in M30's central region. Figure ^ is a composite of the {B — V, V) CMDs in each 
of the eight radial bins. The morphology of the stellar distribution in this CMD is qualitatively 
similar to that in the {B - V, B) CMD (Fig. g), except that both the RGB and HB are slightly 
more extended and vertical than in the latter CMD. There are no striking differences in the mix 
of stellar types from one annulus to another, aside from the obvious excess of BSSs in the inner 



bins (see also Sec. \3.2.^ and Table y). There is an indication that the central regions are deficient 
in bright RGB/AGB stars (V ^ 16); the innermost bin also appears to be slightly deficient in 
HB stars (the degree and significance of these effects are quantified below). The scatter due to 
photometric error decreases, and the degree of completeness for stars fainter than about V = 20 
increases, as one moves from the dense central region to the less crowded cluster outskirts. 



Table 4. Relative radial distribution of evolved stellar populations in M30. 





Radial Bins 


a 


Stellar Subsample 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


n lb lo" 


Post MSTO stars {V < 18.6) 


160 


159 


155 


159 


157 


157 


153 


160 


158 ± 12.5 


Bright stars {V < Vub + 2) 


36 


43 


40 


43 


47 


40 


36 


47 


41.5 ±6.44 


All stars {V ^ 23) 


645 


915 


939 


1099 


1278 


1449 


1647 


1535 


1188 ± 34.5 


Faint RGB/Subgiants 


93 


86 


87 


85 


87 


104 


96 


85 


90 ±9.51 


Bright RGB/AGB 


5 


8 


9 


16 


16 


13 


15 


16 


12.3 ±3.50 


HB 


8 


16 


13 


12 


13 


11 


9 


16 


12.3 ±3.50 


BSS 


21 


10 


3 


3 


5 





3 


3 


6.0 ±2.45 



^ (1) r < 5'.'00, (2) 5'.'00 <r< 9'.'80, (3) 9'.'80 <r< 15'.'41, (4) 15'.'41 < r < 23'.'2, 

(5) 23'.'2 < r < 35'.'8, (6) 35'.'8 < r < 51", (7) 51" < r < 71", (8) 71" < r < 130". 
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Fig. 8. — A composite oi V vs B — V color-magnitude diagrams in eight radial bins around the 
center of M30. The limiting radii have been chosen to ensure that the bins contain equal numbers 
of stars with V < 18.6. The outer five bins (7^4-7^8) represent partial annuli. The innermost bin 
contains an excess of blue stragglers and appears to be deficient in bright RGB stars relative to 
the bins farther out. The photometric accuracy improves and degree of completeness for faint stars 
increases with increasing radius. 
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Table ^ lists the number of stars of each type in the eight radial bins. The division of stars 
into various stellar types is defined by the dashed lines in Fig. |^. The total number of stars in each 
radial bin increases towards larger radii due to the increased probability of detecting faint MS 
stars {V ^ 20) in the relatively sparse outer parts. Faint RGB stars and subgiants make up the 
bulk of the complete reference sample used to define the limiting radii, so it is not surprising that 
their numbers are roughly constant across the bins. There is a 2.1a depletion of bright RGB/AGB 
stars in the inner bin compared to the average of all bins. The inner depletion is slightly more 
significant (2.6a) relative to the average of the outer four bins (15.0 it 3.9) rather than the overall 
average; the overall average is biased low due to the RGB depletion in the center. The bright RGB 
sample size is too small to tell whether the central deficiency is stronger for the brightest giants 
(those with V < Vhb or F ^ 15) than for the full set of bright RGB stars (V ^ 16 as defined in 
Sec. ^ ); the effect is less significant for the V ^ 15 subset than for the full set, but this is due 
to the increased Poisson error associated with the brighter subset. The slight apparent depletion 
of HB stars in the central 5" compared to the overall average is not statistically significant (1.2a 
effect). 

Figure ^ shows the cumulative radial distributions of different types of stars in the core of 
M30 over the full area of the WFPC2 mosaic (r < 130"). Figure ^ shows the cumulative radial 
distribution for the inner r < 25" region of the cluster. A two-sided Kolmogorov-Smirnov (KS) 
test is used to measure the probability that a given pair of stellar types is drawn from the 
same underlying radial distribution function. We use the faint RGB sample as the comparison 
population for KS tests as this phase of stellar evolution is relatively uncomplicated and well 
understood. Faint RGB, bright RGB/AGB, and HB stars are expected to have the same radial 
distribution because: (1) stars evolving "normally" should proceed from faint RGB/subgiant to 
bright RGB/AGB to HB without a significant change in mass, and (2) even if the masses of 
these different stellar types were slightly different because of mass loss during the bright RGB 
phase, mass segregation proceeds on a dynamical timescale which is longer than the evolutionary 
timescale in these phases. 

The KS probability that the BSSs are drawn from the same radial distribution as the 
faint RGB/subgiant population is less than 10^''. The cumulative radial distribution plots show 
that the BSSs are strongly centrally concentrated with respect to all the other stellar types, and 
that this trend is strong even within the inner 25" of the cluster (Figs. |9| and 10), in keeping with 



the specific frequency vs radius calculation in Sec. and Table |. The bright RGB/AGB 



and faint RGB/subgiant populations have different radial distributions, with only a 5.4% KS 
probability of being drawn from the same distribution. In contrast, the radial distributions of the 
HB and faint RGB/subgiant stars are consistent with each other. 

Previous studies of M30 have noted that its central region is deficient in bright RGB stars 



relative to its outer region (cf. Piotto et al. 1988|). The cause of the bright RGB depletion in 



M30 and other post-core-collapse clusters remains something of a mystery. Djorgovski et al. 
(1991) examine a variety of physical mechanisms — tidal interactions, stellar collisions, Roche 
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Fig. 9. — The cumulative radial distribution of various types of post main sequence stars in the 
WFPC2 image of M30 (r < 130"). Only every tenth faint RGB/subgiant star is shown for the 
sake of clarity. The numbers within parentheses indicate the number of stars of each type. The 
BSSs arc significantly more centrally concentrated and the bright RGB/AGB stars appear to be 
somewhat less centrally concentrated than the faint RGB/subgiant normalizing population. 



- 28 - 




5 10 15 20 

Radius: r [arcsec] 

Fig. 10. — Same as Fig. ^ for the inner r < 23" of M30. Only every fifth faint RGB/subgiant star 
is shown for clarity. The trends observed in the full sample (Fig. |9|) persist in the inner region 
of M30 — both the strong central concentration of BSSs as well as the slight central deficiency of 
bright RGB/AGB stars. 
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lobe overflow of an evolving red giant in a tight binary system, etc. — but find that none of these 
provides a convincing explanation of the observations. The general consensus appears to be that 
binaries play some role in causing the central depletion of bright RGB stars; this has motivated 
several detailed, but as yet inconclusive, theoretical investigations of the effect of binaries on post 
MS evolution (cf. [Taam fc Lin 1992| ; P'Antona et al. 1995| ). 



3.3.2. Which Stellar Populations are Responsible for the Radial Color Gradient Near M30's 

Center? 

In order to study the radial color gradient in the integrated starlight of the cluster, we 
construct mosaic images from the WFPC2 data in the F439W {B) and F555W {V) bands (greyscale 
image in Fig. |l]) using the metric task in iraf/stsdas which corrects for distortion/aberration 
in each CCD image and inter-CCD rotations and translations. Standard aperture photometry 
routines (with sky subtraction) are used to measure the integrated flux in each band in each of 
the eight radial bins defined in Sec. The background "sky" level is estimated from relatively 
sparse areas of the mosaic image. Sky counts make up as much as 22% of the total counts in 
the outermost radial bin (71" < r < 130") of the F439W image; the fractional contribution of 
sky counts is lower in the F555W band, and decreases sharply at smaller radii in both bands. 
The integrated cluster fluxes are converted to B- and y-band surface brightnesses using the 
photometric calibration parameters (zeropoint and color transformation) described in Sec. |2.2| . 
The "effective" radius of each bin is taken to be the median radial distance from the cluster 
center of all evolved stars in that bin. This is more representative than the simple average or 
area-weighted average of the limiting radii of the bin; the distinction is especially important for 
the innermost bin which covers the largest dynamic range in radius, and for the outermost bin 
whose outer edge is defined by the outline of the WFPC2 mosaic image and this causes the actual 
number of stars detected to be a strongly decreasing function of radius. 

Figure ^ shows the integrated B — V color of M30 as a function of radius (bold line with 
dots in both panels). Our measurements are in good agreement with previous measurements of 
the color gradient of the overall cluster light ( piotto et al. 1988 ; Burgarella Sz Buat 1996| ) and may 



be summarized as: d{B — V)/dlog{r) = 0.3 mag/dex with a color of {B — V)ref ~ 0.59 at a radius 
of rref = 30". 

We next investigate which of M30's stellar types is/are responsible for the observed color 
gradient. If the color gradient were caused by a single population of stars, the gradient would 
disappear when light from that population in each band is redistributed "uniformly" throughout 
the cluster. We use the radial distribution of the integrated B + V stellar flux to define what 
is meant by "uniform" redistribution. Before summing the fluxes in the two bands, the relative 
flux normalization is chosen such that faint RGB stars make an equal contribution in B and V; 
these stars have the same color as the mean color of the integrated cluster starlight over the area 
of the WFPC2 image: {B - V)m30 ~ 0.65. For each population (BSS, HB, bright RGB/AGB, 
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Radius: r [arcsec] 



Fig. 11. — The integrated B — V color of M30's starlight as a function of radius (bold solid line 
with dots) . Each of the other curves illustrates the effect of "uniformly" redistributing the light of 
a particular stellar type (see Sec. 3. 3.^ . The mean B — V color of M30 over the area of the WFPC2 
image (r ^ 2') is indicated in the lower panel. The integrated cluster starlight is about 0.4 mag 
bluer in i? — y in the inner part than in the outer part. Bright RGB/AGB stars are responsible 
for about 0.15 mag of this color gradient, and HB stars for about 0.05 mag. The evolved stars 
taken together are responsible for about half of the overall color gradient, while main sequence 
stars {B > 18.75) are responsible for the remainder of the gradient. 
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faint RGB) in turn, we subtract the actual flux of that population in every radial bin, and then 
redistribute the total amount of flux removed in a uniform and identical fashion in both bands. 
The step of adding the light back in is critical as it dilutes any remaining color gradient by the 
appropriate amount, thereby allowing us to quantify the contribution of a given stellar type. 

The effect of redistributing the light of a single evolved stellar population on the overall 



color gradient is shown by the various curves in the upper panel of Fig. 11. It is clear that 
bright RGB/AGB stars affect the color gradient more than any other evolved stellar type, but 
account for only about a third of the observed gradient, or only d{B — V)/d\og{r) = 0.1 mag/dex. 
The observed lack of smoothness in M30's integrated B — V profile is a reflection of the stochastic 
nature of the bright RGB/AGB flux component: a relatively small number of stars contributes 
a significant fraction of the total light. Redistribution of the light of bright RGB/AGB stars 



results in a considerably smoother color gradient (short dashed line in upper panel of Fig. 11). 
The population of blue HB stars, long a popular suspect for the central bluing trend, do not affect 
the the overall color gradient of the cluster light very much: d{B — V) / d\og{r) = 0.05 mag/dex. 
Another popular suspect, the BSSs only affect the color by 0.02 mag in the innermost bin. 

Although quite numerous near the center, the BSS sample is not a significant contributor 
to the overall gradient. A simple calculation verifies this: bright RGB stars are about 3 times 
as numerous as BSSs in the center of M30 (see -Frss in Sec. |3. 2. ^ ) ; they have mean colors of 
{B — y)bright RGB ~ 0.6 and {B — y)BSS ~ 0, a difference of A(i? — V) = 0.6 mag; the typical BSS 
is about 2 mag fainter than a bright RGB star (the ratio of their fluxes is ~ 15%); bright RGB 



stars contribute about 50% of the total cluster light (see Sec. 3^); the fractional BSS portion of the 
integrated stellar flux in the central bin is 0.33 x 0.15 x 0.5 = 2.5%; thus, the concentration of BSSs 
in the innermost bin is expected to affect the overall B — V color by only 0.025 x 0.6 = 0.015 mag. 

The lower panel of Fig. |ll| shows that a smooth color gradient of d{B — V)/d\og{r) = 
0.15 mag/dex remains when the flux from all evolved stars is redistributed. This remaining color 
gradient must be due to the unevolved stars (with B > 18.75) in M30. We estimate the combined 
flux of all cluster MS stars — both resolved stars with V <^ 21 as well as those that are fainter 
than the detection threshold — by subtracting the contribution of all evolved stars from the total 
cluster flux in each bin. By redistributing the MS light uniformly, we verify that this population 
does indeed account for approximately half of the overall color gradient observed (dashed line in 



the lower panel of Fig. 11). The contribution of MS stars to M30's radial color gradient is best 



explained in terms of mass segregation ( Pryor et al. 1986 ; Bolte 1989| ). There is direct evidence 



that the fractional flux of MS stars increases outwards across the eight radial bins (see Sec. "iA 
below). Moreover, the mean color of the MS light is expected to become redder with increasing 
radius because of an increasing fraction of lower MS stars whose colors are redder than the 
cluster's mean B — V color of 0.65. These faint red stars are mostly unresolved in the WFPC2 
images; the resolved MS stars tend to be the brighter ones just below the MSTO and these have 
bluer colors (Fig. 0). 
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Djorgovski & Piotto (1993) suggested that M30's color gradient is caused by depletion of giants 
in center of the cluster, while Burgarella & Buat (1996) concluded that the color gradient cannot 
be explained by the evolved and resolved stellar populations alone. We find that bright RGB/AGB 
stars are responsible for only about a third of the overall color gradient and that about half of 
the gradient is caused by stars below the MSTO. There is a 5% chance that the central depletion 
of bright RGB stars, and the associated color gradient of d{B — V)/dlog(r) = 0.1 mag/dex, is a 
statistical fluke (Sec. 3.3.1 ). The contribution of HB stars to the gradient is more likely to be a 
result of Poisson fluctuations. The color gradient due to the mass segregation of MS stars, about 
half the overall gradient, is highly significant in a statistical sense. 



3.4. The Contribution of Various Stellar Types to the Integrated Starlight 



Figure 12 shows the V surface brightness as a function of radius for various stellar populations 
as well as for the total cluster light (bold line with dots). Bright RGB/AGB (defined by < 16.6 
for this study) and HB stars taken together account for most of the light of evolved stars, especially 
beyond r > 20". The central depletion of bright RGB/AGB stars relative to faint RGB stars 
is apparent: the integrated bright RGB/AGB/HB flux is 4 times that of the faint RGB flux in 
the outer parts, whereas their fluxes are comparable in the innermost bin. The MS component 
(total flux minus integrated flux of all evolved stars) is about a third of the total flux in the inner 
two bins and increases to about half the total flux in the outermost bin, likely a result of mass 
segregation. The "true" degree of segregation of the MS component is slightly larger than this 
observed fractional increase, since the effect is diluted by the central depletion of bright RGB/AGB 
stars. For example, the ratio of integrated MS flux to faint RGB flux increases from 1.5-2 in the 
inner part (r ^ 10") to about 4-5 beyond r > 50". 

The surface brightness profile of M30's integrated starlight displays significant curvature in a 
log-log plot (Fig. [l^). The slope interior to r = 10" is about a = —0.5, where a is the power law 
index of the projected density profile: iy[r) oc r"^. The faint RGB component's brightness profile 
slope in this same inner region is a ~ —1; this measure is unaffected by the central depletion of 
bright RGB stars and is in good agreement with the number-weighted estimate of M30's density 



profile (Paper IV). The profile slope for the total light as well as for each stellar type steepens to 
about a = —2 beyond r > 40", the asymptotic power law slope of a King profile for r Vcore 
(iKing 19621). 



3.5. Luminosity Function 

The observed distribution of red giant luminosities in a globular cluster, while relatively 
easy to measure, can serve as a direct and powerful diagnostic of the rate of fuel consumption 
as a function of the evolutionary stage of the star from the beginning of the post-main-sequence 
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Fig. 12. — The V band surface brightness as a function of radius in M30 for the integrated 
cluster starhght (bold solid line with dots), as well the separate contributions of various stellar 
types: all post main sequence stars, bright RGB/AGB/HB stars, and faint RGB/BSSs. The main 
sequence component (total minus evolved stars) increases radially outwards relative to the total 
light, and especially in relation to the faint RGB component, as a result of mass segregation. 
Bright RGB/AGB/HB stars contribute 3-4 times as much flux as faint RGB stars, except in the 
innermost bin which is deficient in bright RGB stars. The profile of the integrated light has a slope 
consistent with an a = —0.5 power law for r < 10" (a = — 1 for the faint RGB light) and steepens 
to a = -2 for r > 40". 



- 34 - 



phase (development of hydrogen burning sheh) to the tip of the RGB phase (hehum flash). The 
luminosity of an RGB star is expected to be determined principally by the mass of its helium core 
( Refsdal fc Weigart 1970| ). Detailed stellar evolutionary models do indeed confirm that the gross 



shape of the LF of red giants (those significantly brighter than the subgiant branch) is largely 
independent of the metallicity and age of the cluster (Bergbusch & VandenBerg 1992). 



Figure compares the observed M30 F-band stellar LF to model isochrones computed 
by Bergbusch & VandenBerg (1992) with abundances of [Fe/H] = -2.03, [0/Fe]= +0.70, 
and Y = 0.235. The adopted composition is consistent with abundance estimates for the 



cluster (cf. |Zinn &: West 1984 ; Carretta fc Gratton 1997 ; Sandquist et al. 1998 ). An apparent 



distance modulus of {V — Mv)app = 15.10, based on fitting subdwarfs with Hipparcos distance 
measurements to the cluster MS ( [Reid 19*97 ) and line-of-sight reddening of Eb-v = 0.05 ( |Zinn 



19801 ; IBurstein fc Heiles 1982| ; [Richer et al. 1988| ), is adopted. Most AGB stars (all but the reddest 



ones) and all HB stars are excluded in computing M30's LF since such stars are absent from 
the Bergbusch & VandenBerg model; the remaining AGB stars make a negligible perturbation 
on the LF, as do BSSs (^ 10% effect for the bins around My = +2 to +3). No completeness 
correction has been applied; incompleteness is unimportant across most of the brightness range 
shown in Fig. |l^ {My < +5.5), as the stars fainter than the MSTO (solid dots) are drawn only 
from the outer two annuli (51" < r < 130") where crowding is not too severe. The relatively rare 
bright stars comprising the My < +2.5 portion of the LF are counted over the full WFPC2 image 
(r < 130"); the bright and faint parts of the LF are matched to each other using counts of all stars 
with V < Fmsto (My < +3.5). 



The Bergbusch & VandenBerg model LFs shown in Fig. 13 are normalized to match the 
data over the range 19.5 ^ V ^ 20.0 (+4.5 ^ My ^ +5.0). The 10 Gyr and 8 Gyr isochrones 
provide acceptable fits to M30's MS stars in the range +3.5 < My < +5.5 (incompleteness 
sets in for My ^ +6 or ^ 21) for initial mass function slopes in the range x = to x = 1, 
where x is defined by the usual relation: dN{M) oc M~^^~^^^ dM . The portion of the LF with 
My < +5.5 corresponds to a very narrow range in stellar mass (0.73 Mq < M < 0.90 for a 
10 Gyr isochrone) so that the shape of this portion of the model LF is almost independent of the 



mass function slope x. Based on studies extending further down the MS beyond r > 2' ( Bolte 
19941 ; gandquist et al. 1998|) , a mass function slope of x = 1 is quite reasonable for the r ~ 1' - 2' 
region of the cluster from which the faint stars in Fig. ^ are drawn. While the isochrones are 
normalized to the MS portion of M30's LF, the 10 Gyr isochrone also fits the subgiant portion 
well. The MSTO luminosity for the 8 Gyr isochrone is noticeably brighter than in M30, and this is 
manifested as an excess in the predicted number of subgiants across the range +2.5 < My < +3.5. 

The number of RGB stars with My < +2 (y < 17.1) found in M30 is 220, 30% higher than 
predicted by the appropriately normalized model LFs (170). The deficit is significant at the 3.8cr 
level. The effect is independent of RGB luminosity for My < +2. Earlier ground-based studies 
have noted a similar mismatch between observation and standard stellar evolutionary models in 



the r > 2' region of M30 (|Bolte 1989| ; [Sandquist et al. 1998| ). The "RGB excess" phenomenon 
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Fig. 13. — The F-band stellar luminosity function of M30 (bold dots with Poisson error bars and 
bin width markers) compared to the Bergbuseh &: VandenBerg (1992) models for 8 Gyr and 10 Gyr 
metal-poor populations for various mass function slopes x (lines) . The models are normahzed to the 
observed number of stars in the range My = +4.5 to +5.0. The data points for bright {My < +2.5) 
stars are based on the entire WFPC2 mosaic; the faint star data (bold dots) are based only on 
the relatively sparse outer radial bins (r ^ 50"). The turnover in the counts for My > 6 is due to 
incompleteness; the effect of incompleteness is stronger and sets in at brighter magnitudes closer 
to the cluster center (open circles). The models underpredict the number of bright red giants by 
30%. 
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is not unique to M30 — the metal poor clusters M5 (Sandquist et al. 1996) and M13 (Paper VI) 
display an RGB excess relative to theoretical models, but unlike the case of M30, the excess in 
these two clusters is restricted only to giants brighter than the RGB bump or HB level {My ^ 0). 

One possible way of resolving the discrepancy between M30's observed RGB LF and 
theoretical LFs is by invoking the effect of internal core rotation in stellar models ( Pearson et al. 



1995; VandenBerg et al. 1998). This tends to decrease the rate of RGB evolution (by expanding 



the core and lowering the shell temperature) and to thus increase the number of giants in any 
luminosity bin at a given instant. Core rotation models have only been constructed, however, for 
giants fainter than My = 0; thus, these models do not directly address the RGB excesses observed 
in M5 and M13. A more comprehensive solution based on "deep mixing" has recently been 



explored by Langer et al. (1998) . Deep mixing has already been invoked to explain: (1) chemical 
abundance anomalies in some bright cluster giants in terms of processing of envelope material in 
the hotter energy generating shell (cf. [Kraft 199^ ); and (2) anomalies in the HB color (second 
parameter problem) and brightness in terms of introduction of helium rich material into the 
envelope ( ^weigart 1997| ). Langer et al. propose that the "extra" fuel mixed into the hydrogen 
burning shell might also account for the overabundance of RGB stars seen in certain clusters. 

The observed excess of red giants relative to the standard model is not an artifact resulting 
from incompleteness in the faint star counts. As discussed above, the shape of the MS LF for 
+3.5 < My < +5.5 is a good match to model LFs, so any incompleteness in the star counts would 
have to be independent of apparent magnitude over this range (18.5 < V < 20.5). Our image 
simulations indicate that the sample is complete down to y ~ 21 over the r ~ 1' - 2' region from 
which this part of the LF is constructed (Sec. ^^). Moreover, incompleteness tends to be a strong 
function of apparent magnitude; this is clearly seen in the crowded inner region (5" < r < 15") of 
M30 where incompleteness sets in at brighter magnitudes (open circles in Fig. |l^; see also Fig. |2|). 

The conclusion about M30's RGB excess is essentially independent of the adopted cluster 
distance modulus (m — M) and choice of isochrone age t (Fig. |l^. For any reasonable value of 



(m — M), t can be chosen so that the model LF (Bergbusch &: VandenBerg 1992) fits the MS and 



subgiant "bump" portions of M30's LF, but the model underpredicts the number of RGB stars. 
Even if the assigned value of t is too large (too small) for a given distance modulus, the result is 
an apparent excess (deficiency) in M30's subgiant population relative to the model, provided the 
model LF is normalized to the MS portion of M30's LF; the degree of M30's RGB excess though 
is unaffected by such a mismatch between t and (m — M). In spite of recent improvements in 
the measurement of subdwarf parallaxes with Hipparcos, the distance modulus of M30 remains 
uncertain at the level of A(m — M) ^ 0.2 mag, corresponding to an uncertainty in the cluster 
age of At ^ 2 Gyr ( |Reid 1997] ; [Gratton et al. 1997] ; [Font et al. 199^ ; pandquist et al. 1998| ). The 



distance scale we have adopted, {V — -/Vfy)app = 15.10 ( |Reid 1997 ) happens to be at the high end 



of this range, so that the corresponding best-fit age, t ~ 10 Gyr, is lower than in other studies. 
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4. Conclusions 

1. Hubble Space Telescope Wide Field/Planetary Camera 2 images of the dense globular cluster 
M30 (NGC 7099) in the F555W, F439W, and F336W filters have been analyzed. Accurate 
stellar positions and photometry in the (Johnson) UBV bands are presented for 9507 stars 
(and BV photometry for an additional 433 stars) within a projected distance of 130" from 
the cluster center. Color-magnitude diagrams based on the UBV bands are presented 
showing clearly distinguished sequences of stellar types: red giant branch, (blue) horizontal 
branch, subgiant branch, blue straggler sequence, main sequence turnoff, and main sequence. 
Typical photometric errors for main sequence stars {B = 19-20) range from ~ 0.1 mag for 
the central 10" of M30 to ~ 0.05 mag in the rest of the cluster. Incompleteness sets in a 
little fainter than the main sequence turnoff (V ^ 19) for r < 5" but the sample is complete 
to progressively fainter magnitudes at larger radii {V ~ 20.5 for r ^ 20"). 

2. An unusual star is found 1"2 SSW of the cluster center. Its very blue color {B — V = —0.97, 
B = 18.6) indicates that it might possibly be a cataclysmic variable. 

3. Forty-eight blue straggler candidates are identified in the cluster on the basis oi a {B — V, B) 
color-magnitude diagram; the BSS classification for some of these is tentative. The specific 
frequency of BSS in M30, Fbss = iV(BSS)/A^(y < Vhb + 2, is 0.39 ± 0.08 in the inner 10" 
and 0.14 it 0.02 over the full WFPC2 image (r < 130"). The central BSS specific frequency 
in M30 is the highest measured for any globular cluster to date. The BSSs are strongly 
centrally concentrated relative to other evolved stars (> 99% significance). The observed 
M30 BSS luminosity function is compared to theoretical predictions: although the number 
of BSS is too small to draw definite conclusions, the collisional formation model seems to be 
a better match to the data than the primordial binary merger model. 

4. The abundance of bright red giant branch stars and asymptotic giant branch stars with 
y ^ 16 (relative to fainter RGB stars and subgiants) appear to be a factor of 2-3 lower in 
the central 15" of M30 compared to further out in the cluster (2-2. 5o" effect). Horizontal 
branch stars, the evolutionary descendants of bright RGB stars, show no significant central 
depletion relative to the distribution of faint RGB stars and subgiants. 

5. The B — V color of the integrated cluster light in M30 varies from 0.82 around r ~ l' to 
0.45 within r < 10", a radial color gradient of d{B — V)/dlog(r) = 0.3 mag/dex. The 
central deficiency of bright red giant / asymptotic giant branch stars is responsible for about 
a third of the color gradient; all evolved stellar populations taken together (RGB, AGB, 
HB, BSS, subgiants) are responsible for about half of the overall color gradient. Mass 
segregation of faint red main sequence stars results in a smooth B — V color gradient of 
about 0.15 mag/dex, corresponding to half the observed gradient. 

6. The F-band surface brightness profile of the integrated starlight in M30 has a power law 
slope of a = —0.5 in the inner 10", and the slope is a = — 1 for the light of faint RGB 
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stars alone; the latter is identical to the number-weighted measure of the projected density 



profile slope found in Paper IV. The difference between the central slopes of the overall vs 
faint RGB light is a result of the central depletion of bright RGB stars. The integrated 
surface brightness profile slope steepens to a = —2 for r > 40". Main sequence stars 
contribute nearly 40%, bright RGB/AGB/HB stars with V ^ 16 contribute 30%, and 
faint RGB stars contribute about 30% of the overall cluster brightness in the central 5" of 
M30; the corresponding fractions are 50%, 40%, and 10% at larger radii (r ~ 1'). 

7. The stellar luminosity function in the inner 2' of M30 shows an anomalously high abundance 
of red giants (relative to subgiant, main sequence turnoff, and main sequence stars), as 
has been previously noted further out in the cluster. Giants with My < +2 in M30 are 
overabundant by 30% (3. Sex significance) relative to the prediction of suitably normalized 
standard stellar evolution models. 
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Appendix: Table of Stellar Photometry and Positions 

The complete version of Table |l] contains (Johnson) U-, B-, IZ-band photometry for 9507 stars 
with F ^ 21 in the central 5.1 arcmin^ of M30 (r < 130") covered by the HST WFPC2 image, and 
B and V for an additional 433 stars that were not matched in the U band; a subset containing the 
40 brightest stars within r < 10" of M30's center is presented in the paper. A computer-readable 
version of the complete table is available in the electronic edition of the Journal by link to a 
permanent database. It may also be obtained via anonymous ftp as follows: 

• ftp eku. sns . ias . edu (login as: anonymous) 

• cd pub/GL0BULAR_CLUSTERS/m30-ppr7 

• get m30fulltabl.txt 

or by contacting P.G.. 

Details of the photometric procedure, a description of the astrometric convention, and 



estimates of the completeness and photometric accuracy may be found in Sec. |2.2| , |2.3| , and [2A , 
respectively. We recommend that the astrometry in Table |l| be matched empirically against other 
data sets to correct for possible systematic errors in the scale, translation, and rotation of the 
coordinate system we have adopted. 
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